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We study the correction to the radiative decay width due to finite light quark mass. Using lattice 
QCD with the overlap fermion formulation, we calculate the three-point function of the form 
{PVixVv) in the (Euclidean) momentum space, which corresponds to the tt*^ ^> 7*7* amplitude. To 
fit the lattice data, we use two different modifications of vector meson dominance (VMD) ansatz. 
One is a combined form of VMD with the next-to-leading order (NLO) chiral perturbation theory 
(ChPT), and the other is a resummed form of pion-loop diagrams. We extract one of the low 
energy constants in NLO ChPT, and estimate tt" jj decay width. 
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1. Introduction 

We present an update on our study [P of the 

TjO — )• yy decay width using lattice QCD calcu- 
lation of its off-shell amplitude. In the limit of massless up and down quarks, the corresponding 
amplitude is completely determined by the triangle anomaly which is valid to all orders in 

perturbation theory [||]. Once the quark mass is introduced, the correction term arises, which can be 
classified most conveniently in the effective field theory framework [|5j, ^. The most recent analysis 
based on the next-to-next-to-leading order chiral perturbation theory including the isospin breaking 
effect quoted Y^^^ = 8.09(11) eV, which marginally agrees with the recent experimental data 
Fj^a^yy = 7.82(22) eV [Q]. The experimental error is expected to be reduced by a factor of 2 in 
the near future, and the improvement of the theoretical prediction would become more important. 
Presently, the theoretical estimate is limited by the lack of knowledge of the low energy constants 
in the effective theory, for which the lattice input is most desirable. 

We calculate the ;r*' ^ 7*7* (two photons are off-shell) amplitude in lattice QCD from the 
(PV^Vv)-type three-point function in the momentum space. In our previous report [|]], we showed 
that the tt" ^ 7*7* amplitude calculated on the lattice using the overlap fermion formulation is 
consistent with the expectation from the axial anomaly, which is 1/(471^), after extrapolating to 
the limit of on-shell photons and massless quarks. The overlap fermion is most suitable for this 
study, because it has an exact U(1)a chiral symmetry at the Lagrangian level that is violated by the 
quantum effect. Namely, it has the same symmetry structure as in the continuum theory even at 
finite lattice spacings. 

We work on a two-flavor QCD ensemble of a 16^ x 32 lattice with the lattice cutoff l/c? = 
1.667(17) GeV [p^. Dynamical quarks are also described by the overlap fermion. Topological 
charge Q in this ensemble is fixed at 2 = 0, which is a source of the finite size effect JTl]]. There 
are four different quark masses, rriq = 0.015, 0.025, 0.035, and 0.050 (in the lattice unit), which 
cover the range between ms/6 and ms/2 with physical strange quark mass m^. 

In the analysis presented last year [|]], we used the naive vector meson dominance (VMD) 
ansatz to fit the lattice data and to extrapolate to the Umit of on-shell photons. In the present report, 
we update this work including an analysis using other functional forms motivated by the chiral 
effective theories. We then present our preliminary result for the pion decay width. 

Once the tt" 7* 7* amplitude is determined including its off-shell form factor, it may be used 
to estimate the size of the light-by-light scattering contribution to the muon g — 2, which occurs 
through a (probably dominant) virtual process 77* — )■ 71° — > 7*7*. 



2. n^' ^ Y*y* transition form factor 

The 71° — )• 7*7* transition form factor fji0yt.y*{pi,p2) is defined as a matrix element of two 
electromagnetic (vector) currents V™ between the pion state and the vacuum: 



2 



Two-photon decay of 7Z from two-flavor lattice QCD E. Shintani 



where pi and p2 denotes the off-shell photon momenta and q = —p\ — P2- On the lattice we extract 
fjfiy*y^{pi,P2) from a three-point function 



= 2m,tv[z'Ql] / ^ ^'-'e--'fty2ReTr[5,(0,x)rp(x,y)5,(y,3')7v5,(3;,0)7^] \ , (2.2) 

where (• • •) denotes an ensemble average. Sq{x,y) represents the quark propagator Dov{m)^^{x,y) 
for the overlap-Dirac operator D„^,{m) with the mass term m. To achieve the C>(fl') -improvement of 
the pseudoscalar density operator P™'^ contains the structure rp{x,y) = (1 — Doy{0)/MQ)Y5{x,y). 
"tr" means a trace over flavor indices, and "Tr" represents a trace over color and spinor indices. 
The (Euclidean) four-momenta Pi, P2 and Q{= —Pi — P2) are those of two photons and a pion, 
respectively. Their components are limited to discrete values 2nn^/{L^a), with the lattice spacing 
a and the lattice size ^x,y,z = 16 and L, = 32. In obtaining the second hne of ( |2.2D , we assumed that 
the disconnected quark diagrams are negligible. 

With the overlap fermion formulation, the (flavor non-singlet) pseudo-scalar density opera- 
tor P"'^(;c) = ^(;c)t^75[(1 - D„y(0)/Mo)q]{x) satisfies the axial-Ward-Takahashi (AWT) identity 
2m^P"'^ = with the conserved axial-current A^^^ [12], which we numerically confirmed. 



For the electromagnetic current v^o^em ^ ^^^^^ cun-ent defined with the overlap fermion 

field. The electromagnetic charge is assigned as Qe = diag{2/3e, — 1/3^} for up and down quarks. 

Because we use the sequential source method, we have to invert the overlap-Dirac operator for 
each choice of P2, which is inserted at the point y. We restrict P2 to the cases that it has a non-zero 
component only in the temporal direction, P2 = (0,0,0,P2r), with v = 1. Then we construct the 
three-point function by contracting the trace at the point x with a momentum Q = —Pi — P2 inserted. 
In order to minimize the numerical effort and the storage, we also restrict Pi as 

Pi = (0,Piv,0,0) when ^ = 3 or Pi = (0,0,Pu,0) when n=2, (2.3) 

with all possible values of Piy and Pi^. With these choices, we can pick up the non-zero components 



of PlaP2p£^ivap in(gj). 

We then extract the form factor fjiOy'yiPijPi) as 



pl-t(Pl,P2) = G;:7(P2,e)/ {Y^PlaPipe^vap I = -J^^fn^rAPuPl) (2.4) 

up to the contributions from excited states and disconnected diagrams. Saturation of the pseudo- 
scalar channel by the pion field is quite accurate for small momenta, because the first excited state, 
called ;r(1300) in PDG, is much heavier. We numerically confirmed this by looking at the two- 
point function (PP) in the momentum space. 

Since Pi and P2 are defined in the Euclidean space-time, the on-shell condition is realized only 
at the zero momentum for both Pi and P2. In this limit, the three-point function trivially vanishes 
and no information on the form factor could be extracted. In order to obtain the decay width, 
we therefore have to work with slightly off-shell photons and to extrapolate to the on-shell limit 
assuming some functional form for f-^yty* (P\^P2). 
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Figure 1: Vertices contribution to the tt" — > 7*7* decay. The soHd, double-line and waved lines show pion, 
vector meson, and photon propagator, respectively. 



3. Fitansatz 



For the form factor fj^oyy* (Pi^Pi), a naive choice would be that of the vector dominance model 
(VMD) 



7VMD, 



1 



■Gv {Pi,mv)Gv (Pi, my) 



(3.1) 



+ ml 4n^ 

where Gy {P, my ) = my / (P^ + my ) with the vector meson mass my . In the limit of P\ = P2 = 0, 
reduces to l/An^, which is the expectation of the anomaly valid in the massless limit. Note 
that the above approximation is too naive as the same value 1 / (471^ ) is also obtained even at finite 
quark mass, and so we clearly need modifications. 

Apart from this limit, we expect some corrections. For small quark mass and external mo- 
menta, those can be described by the one-loop ChPT. A formula for the cases of off-shell photons 
is available in [[l3|]. We slightly extend the formula by adding some analytic term and use a function 



plloop+VMD(^2^p2_p2) 



1 + ^ ( - ^AM^) +J{ml,Pf)+J{ml,Pt 



f2 
J 71 

256n^ 



+ 



3 -.3'i + ' 3 

p2 p2 

M^2 

{Pf + ml){P^ + ml) 



Pi 



+ ■ 



^2 



Pl + ml Pl + ml 



y3 



(3.2) 



for -47i:'^F{Pi,P2)/{ml/{Q^ + ml)). (Namely, r(m|,P^,P2) becomes 1 in the massless and soft 
photon limit.) Here, the first hne in the above equation corresponds to the one-loop corrections. 
The tadpole diagram (Figure |l|(a)) and the vector-vector correlator diagram (Figure ^b)) lead to 
the following non-analytic functions 



J{ml,P') 



ml ml 



16712 



64712 



m 1 
-1+ln— ^ + 2atanh"' — 



p2 + 6m| 
967r2 



(3.3) 
(3.4) 



respectively, with a = ^J\+4n^J¥^ . The scale pL is the renormalization point that defines the 
counter terms. We set /i = 775 MeV, the physical p meson mass. 
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Figure 2: Bubble diagram of pion loops. Symbols are the same as in Figure |l} Filled circle shows the 
tree-level vertex for the n-n scattering diagram. 



The second and third lines of ( |3.2| ) define the analytic (or counter) terms. Since the lowest 
non-zero momentum ~ 0.43 GeV^ is already quite large to apply the one-loop ChPT formula, 
we modify the analytic terms motivated by VMD. Namely, the term with a single vector meson 
propagator (the forth term) represents the effective interaction of the form Ti^/p (Figure |](c)), while 
the last term represents n^pco (Figure |l](d)). The parameters y\, y2 and ys are to be determined by 
a numerical fit of the lattice data. 

Another fit form we attempted is that of a resummed pion loops. Again motivated by VMD, 
we assume that the bubble diagrams shown in Figure ^ dominate the higher loop corrections and 
eventually produce the vector meson pole after summed up. Then, an expected functional form in 
the channel corresponding to the rho meson would be 

Jn 

L{ml,Plz2) L{ml,Plz2) 
^ \-2Liml,Pf,Z2) ^ l-2Liml,Plz2y ^ ' ' 

L{ml,P\z2) = QA,(m2) +7(m2 + ^p\^, (3.6) 



where z\ and Z2 are fit parameters. 



4. Fitting results 



Here we describe the fit of the lattice data. 

In order to keep the external momenta as small as possible, we restrict the momentum range 
to use, so that only {aP\Y = 0.154 and 0.039 < {aP2Y < 0.154 are included. These correspond to 



the lowest (aPi)^ and the lowest and second lowest {aP2Y in the momentum assignment (23). The 
data points are plotted in Figure ^ as a function of pion mass squared. The plot is normalized so 
that the on-shell limit corresponds to 1 in the chiral limit. So, we can see that the deviation from 
that limit is quite significant (as large as ~ 50%). 

By fitting the lattice data with the functions r"°°P+^'^'^(m2 ,p2,p2) ^nd P'™'^(m2 ,p2 p2^^ 
we obtain the solid and dashed curves shown in Figure ^ respectively. Here we input the measured 
value of my from two-point function {VV) into r"''°P+^'^°(m|,ff ,ff ). The naive VMD is also 
plotted by thin dashed curves for comparison. In the plot, one can see that both fit curves with the 
modified VMD show a better agreement with the lattice data than the naive VMD does, especially 
for the second lowest momentum. 
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o (aP/=0.03S 

□ (aP/=0.J54 
—- VMD 
— 1-loop+VMD 
--■ Resum 



(aPj) =0.154 




(am ) 

Figure 3: ;r*' — > 7*7* transition form factor as a function of pion mass squared. Different symbol denotes 
the different (0^2)^, and (aPi)^ is fixed in the lowest value, 0.154. Thick solid and dashed curves are results 
in fit function pHoop+VMD pesum respectively. Thin dashed curve shows the naive VMD results as the 
comparison with them. 









3^3 


1-Ioop + VMD 


-0.00030(20) GeV-2 


-0.0054(6) 


1.65(48) 




Zl(= W3) 






Resummation 


-0.00120(12) GeV-2 


-0.0271(37) 





Table 1: Fitting results with the function r^^°°'P+^'^°{ml,Pf,P^) and with F'='™(ot2 ,f f^pf). yi{zi) corre- 
sponds to LECs, W3. 



Fit parameters are listed in Table |I[ A comparison with the phenomenological analysis in 
[13] is possible for the parameter controlling the pion mass dependence W3, which corresponds 
to yi and z\ in our functional form. Unfortunately, we find a significant discrepancy between yi 
and zi- This suggests that the extraction of the low energy constants in the chiral effective theory 
from our calculation suffers from large systematic error. That is due to several possible reasons. 
One of them appears to be due to the large momentum we took for the off-shell photons to apply 
the one-loop ChPT. we estimate this parameter as W3 = — 0.00075(23)(63)GeV^^, where the first 
error is statistical the second error denotes the systematic error coming from the ambiguity of the 



fit function. The phenomenological work used a value W3 = 0.0131 |13], which is an order of 
magnitude larger than our result. 

The tt" YY decay width is defined as 

3 2 

1 nOyy - ^^^3 |i on-shell I ^ ^ 
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with QED fine structure constant = 1/137. Inserting the fitting results after taking = P2 =0, 
which is consistent with CHPT formula in on-shell two photon decay, we obtain 

r^o^^ = 7.98(6)(20)eV, (4.2) 

where the first error is statistical and the second is systematic one. While the systematic uncertainty 
is still rather large in our calculation, the result and its error are compatible with the present world 
average, 7.87(47) eV |^, and also with the recent experimental result 7.82(22) eV [^. 

5. Conclusions 

Lattice calculation of the tt*^ — )- 7*7* decay amplitude is feasible for space-like momenta pi, 
P2 and q. We attempt a fit of the lattice data with the functional forms motivated by the NLO 
ChPT, from which we can in principle extract the low energy constants in the effective theory. 
With our current setup, however, the smallest non-zero momentum is already too large to safely 
apply the NLO ChPT, so that we have to introduce some model functions to extend the ChPT form 
to larger momenta. Our proposals partly motivated by the vector meson dominance ansatz describe 
the lattice data very well, but are not an unique choice. The result for the physical decay width thus 
has large systematic effect, but is already compatible with other estimates in the size of uncertainty. 
An extension to larger lattices and an introduction of the twisted boundary condition will improve 
this situation a lot, which is a subject of future studies. 

Numerical calculations are performed on IBM System Blue Gene Solution and Hitachi SRI 1000 
at High Energy Accelerator Research Organization (KEK) under a support of its Large Scale Sim- 
ulation Program (No. 09/10-09). This work is supported by the Grant-in-Aid of the Japanese Min- 
istry of Education (No. 20105002, 21105508, 21764002, 22011012, 22740183 ). 
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